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ABSTRACT 

Studies of massive black holes (BHs) in nearby galactic centers have revealed a tight correlation 
between BH mass and galactic velocity dispersion. In this paper we investigate how the BH mass 
versus velocity dispersion relation and the nuclear luminosity versus velocity dispersion relation 
in QSOs/active galactic nuclei (AGNs) are connected with the BH mass versus velocity dispersion 
relation in local galaxies, through the nuclear luminosity evolution of individual QSOs/ AGNs and 
the mass growth of individual BHs. In the study we ignore the effects of BH mergers and assume 
that the velocity dispersion does not change significantly during and after the nuclear activity 
phase. Using the observed correlation in local galaxies and an assumed form of the QSO/AGN 
luminosity evolution and BH growth, we obtain the simulated observational appearance of the BH 
mass versus velocity dispersion relation in QSOs/ AGNs. The simulation results illustrate how the 
BH accretion history (e.g., the lifetime of nuclear activity and the possibility that QSOs/ AGNs 
accrete at a super-Eddington accretion rate at the early evolutionary stage) can be inferred from 
the difference between the relation in QSOs/ AGNs and that in local galaxies. We also show how 
the difference may be weakened by the flux limit of telescopes. We expect that a large complete 
sample of QSOs/ AGNs with accurate BH mass and velocity dispersion measurements will help 
to quantitatively constrain QSO/AGN luminosity evolution and BH growth models. 

Subject headings: black hole physics - galaxies: active - galaxies: evolution - galaxies: nuclei - 
quasars: general 



1. Introduction 

The existence of massive black holes (BHs) in nearby galactic centers, as a prediction of the widely accepted 
QSO model that QSOs are powered by gas accretion onto massive BHs (e.g., Lynden-Bell 1969; Soltan 
1982; Rees 1984), is now believed to be confirmed (e.g., Kormendy & Richstone 1995; Magorrian et al. 1998; 
Gebhardt et al. 2003; Pinkney et al. 2003). Studies of the BHs in nearby galaxies have revealed that BH mass 
in nearby galactic centers is tightly correlated with galactic velocity dispersion (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000a; Tremaine et al. 2002) and is also (less tightly) correlated with the luminosity (or mass) 
of elliptical galaxies or bulges of spiral/SO galaxies (e.g., Kormendy & Gebhardt 2001). These correlations 
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suggest a close link between the formation and evolution of BHs and their host galaxies. However, why these 
correlations exist and whether they also exist in distant galaxies have not yet had definite answers. 

The purpose of this paper is to investigate the relationship between BH mass and velocity dispersion in 
QSOs/active galactic nuclei (AGNs). It is important to study this in QSOs/AGNs at least for the following 
two reasons: 

1. Various BH growth models have been proposed to explain the origin of the tight correlation between 
BH mass and velocity dispersion in nearby galaxies (e.g., Silk & Rees 1998; Fabian 1999; Blandford 
1999; Ostriker 2000; Haehnelt & Kauffmann 2000; Burkert & Silk 2001; Adams et al. 2001; Islam 
et al. 2003; King 2003), including accretion of either baryonic gas or non-baryonic dark matter onto 
seed BHs or hierarchical mergers of intermediate-mass BHs (which might be end products of the first 
generation of stars, or Population HI stars) with masses of typically a few hundred M© (e.g., Schneider 
et al. 2002). How can these BH growth models be tested by observations? Among these BH growth 
processes, currently BH mergers are unlikely to be observed directly although detecting gravitational 
wave signals emitted from BH mergers might be possible in the future. The possible observational 
features of accretion of non-baryonic material are still unclear; only accretion of baryonic gas, which 
appears as QSO/AGN phenomena, is detectable and has been extensively studied. We expect that 
the investigation of the BH mass versus velocity dispersion relation in QSOs/AGNs will provide an 
observationally achievable way to understand BH growth and the origin of the correlation of local 
BHs with their host galaxies. As a matter of fact, some observational investigations on the relation in 
QSOs/AGNs have been performed in the past several years (e.g., Laor 1998; Wandel 1999; Gebhardt 
et al. 2000b; Ferrarese et al. 2001; Shields et al. 2003). 

2. Furthermore, according to current observations, the total local BH mass density is consistent with the 
total mass density accreted onto BHs during QSO/AGN phases, which suggests that BH mass growth 
comes mainly from gas accretion during QSO/AGN phases, rather than from accretion of non-baryonic 
material or mergers of intermediate-mass BHs (Yu & Tremaine 2002; Aller & Richstone 2002; Fabian 
2004). Since QSOs/AGNs represent the population of galaxies housing growing BHs, the study of the 
relationship between BH mass and velocity dispersion in QSOs/AGNs can provide valuable information 
on the BH growth history, the evolution of galaxies, and the origin of the tight correlation between BH 
mass and velocity dispersion in local galaxies. 

In this paper we show how the BH mass (or nuclear luminosity) versus velocity dispersion relation in 
QSOs/AGNs is connected with the BH mass versus velocity dispersion relation in local galaxies through 
the nuclear luminosity evolution of individual QSOs/AGNs and mass growth of individual BHs. The basic 
model assumptions are as follows: 

1. BH mergers are assumed to not be important for BH growth. In principle, the mass growth of a BH 
may come; from both gas accretion due to QSO/AGN phases and mergers with other BHs. However, 
currently, the BH merger process and rate are very uncertain. In addition, comparison of the mass 
density distribution in nearby galaxies with that accreted due to QSO phases has shown that BH 
mergers are not necessarily required at least for growth of high-mass (> 10^ M©) BHs (Yu & Tremaine 
2002). 

2. The velocity dispersion of host galaxies is assumed to not change significantly during and after the 
nuclear activity phase. In principle, both nuclear luminosity evolution/BH mass growth and the 
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evolution of host galaxy velocity dispersions contribute to the evolution of the nuclear luminosity /BH 
mass versus velocity dispersion relation. However, ignoring the velocity dispersion evolution may 
greatly simplify our analysis and highlight the effects of nuclear luminosity evolution/BH mass growth. 
Comparison of observations with the predictions obtained by ignoring velocity dispersion evolution 
might also provide clues on the effect of velocity dispersion evolution. In addition, this assumption 
has also been adopted in many other physical models to explain the correlation between BH mass and 
velocity dispersion in nearby galaxies (e.g., Silk & Rees 1998; Fabian 1999; Blandford 1999; Burkert & 
Silk 2001; Adams et al. 2001; King 2003). 

The model includes the following three basic components: the BH mass distribution of local galaxies at a 
given host galaxy velocity dispersion, the curve of the nuclear luminosity evolution/mass growth of individual 
BHs, and the nuclear luminosity/BH mass distribution in QSOs/AGNs at the given velocity dispersion. In 
principle, given any two of them, the last one can be constrained by the model. In this paper we use the 
first two components as model inputs and the third one as the model output. 

The paper is organized as follows: In § 2 we review current observational results on the BH mass 
versus velocity dispersion relation in local galaxies. In § 3 we formulate the connection of the nuclear 
luminosity/BH mass versus velocity dispersion relation in QSOs/AGNs with the BH mass versus velocity 
dispersion relation in nearby galaxies through the nuclear luminosity evolution/mass growth of individual 
BHs. Detailed models on the nuclear luminosity evolution/mass growth of individual BHs are presented in 
§ 4. Applying the observational results of the local BH demography (in § 2) and the assumed QSO luminosity 
evolution and BH growth curves (in § 4) to the analysis made in § 3, we predict the BH mass versus galactic 
velocity dispersion relation in QSOs/AGNs and show the simulation results of its observational appearance 
in § 5. Discussion on current observational results of the relation in QSOs/AGNs is given in § 6. Finally, 
conclusions are summarized in § 7. 

In § 3 readers who are primarily interested in the results of the formulation might want to skip the 
detailed mathematical manipulations and move straight to the model predictions given by equations (11), 
(24), and (26) (with the aid of Table 1). The QSO luminosity evolution or BH growth is incorporated into 
those equations through the lifetime of the nuclear activity rnfe and a distribution function of the probability 
that the progenitor of a local BH with mass Mbh.o (where the subscript "0" represents the current cosmic 
time to) had a mass Mbh or nuclear luminosity L in its nuclear activity history [i.e., P(Mbh|A/bh,o) or 
P(L|Mbh,o) in Table 1]. The observational counterparts of the model predictions are given by equations 
(16) and (30). Comparison of the observations with the model predictions may help to strengthen the 
existing constraints or provide new information on the QSO/AGN luminosity evolution and BH growth. 

In this paper we set the Hubble constant to Hq = lOO/ikms"^, and the cosmological model used is 
{QM,^A,h)= (0.3,0.7,0.65). 

2. The BH mass versus velocity dispersion relation in nearby normal galaxies 

In nearby galaxies BH mass (Mbh.o) is tightly correlated with galactic velocity dispersion (a) (Ferrarese 
& Merritt 2000; Gebhardt et al. 2000a; Tremaine et al. 2002). The logarithm of the BH mass at a given 
velocity dispersion a has a mean value given by 



(log(MBH,o)k) = ^ + 7log((T/200kms-i) 



(1) 
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(Tremaine et al. 2002), where Mbh.o is in units of M©, 7 4.02 ± 0.32, A = 8.18 ± 0.06 has been adjusted 
to our assumed Hubble constant h = 0.65 (see section 2.2 in Yu & Tremaine 2002), and a is the luminosity- 
weighted Une-of-sight velocity dispersion within the effective radius of galaxies. Note that relation (1) is 
fitted in log MBH,o-log cr space. We assume that the distribution in log Mbh,o at a given a is Gaussian, with 
intrinsic standard deviation AiogMeH.o (which is assumed to be independent of a here), and can be written 
as 

(logMBH.o - (logMBH,ok))^ l 

2A2 ■ ^ ' 

logMsH.o 

According to Tremaine et al. (2002), the intrinsic scatter in logMBH.o should not be larger than 0.25-0.3 dex. 
In this paper we set AiogMBH.o = 0.27 dex. 

According to Bayes's theorem, the velocity dispersion probability distribution function (PDF) of nearby 
galaxies with BH mass Mbh,0) p(c|-Wbh,Oj ^0)7 is related to 
p(MBH,o|o-,to) = [ln(10)MBH,o]"^p(logMBH,o|o-,io) by the equation 

,^ n^(a, to)p(AfBH,o|g, ^0) 

p(cr|MBH,0,to) = 7T7 T\ ' 

where riMsH (-^bh,o> ^o) is the local BH mass function (BHMF) defined so that nMBnl-^BH.O) io)dMBH,o 
represents the number density of local BHs with mass in the range Mbh.o ~* -^bh.o + dMBH.o ) {(^^ to ) is the 
velocity dispersion distribution function of local galaxies with massive BHs, and we have nMBH{^BH,o,to) = 
Jna{a, to)p{MBHflW, to)da. 

Besides p(Mbh.o lo", to) above, several other PDFs are also defined below. For comparison of their physical 

mcianings, wv list all of them in Table 1. 



p(logMBH,o|o-,to) = 



27rAi, 



■ exp 



og Mbh.o 



symbol 


physical meaning 


reference 


p(Mbh,o|o-, to) 


PDF of BH mass in local galaxies with velocity dispersion cr 


eq. 2 


p((j|MBH,0,io) 


PDF of velocity dispersion in local galaxies with central BH mass Mbh.o 


eq. 3 


piL\a,t) 


PDF of nuclear luminosity in QSOs/AGNs with velocity dispersion a at 


eqs. 4, 15 




cosmic time t 




p{MBK\(T,t) 


PDF of BH mass in QSOs/AGNs with velocity dispersion a at cosmic 


eqs. 22, 29 




time t 




V{L\a) 


PDF of nuclear luminosity in QSOs/AGNs at all redshifts and with 


eqs. 5, 16 




velocity dispersion a 




V{a\L) 


PDF of velocity dispersion in QSOs/AGNs at all redshifts and with 


eqs. 17, 19 




nuclear luminosity L 






PDF of BH mass in QSOs/AGNs at all redshifts and with velocity dis- 


eqs. 23, 30 




persion (T 




P(L|Mbh,o) 


PDF of nuclear luminosity of the progenitor of a local BH during its 


eq. 8 




nuclear activity phases 




P(Mbh|Mbh,o) 


PDF of BH mass of the progenitor of a local BH during its nuclear 


eq. 25 




activity phases 





Table 1: List of PDFs defined in this paper. 
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3. The demography of QSOs/AGNs 

In this section we analytically link the nuclear luminosity/BH mass versus velocity dispersion relation in 
QSOs/AGNs to the QSO/AGN luminosity evolution, BH growth, and demography of local BHs. We inves- 
tigate the nuclear luminosity versus velocity dispersion relation in QSOs/AGNs in § 3.1 and the BH mass 
versus velocity dispersion relation in QSOs/AGNs in § 3.2. For each relation, we first present the model 
prediction and then show its observational counterpart. 



3.1. The nuclear luminosity versus velocity dispersion relation in QSOs/AGNs 

3.1.1. Model prediction from local BHs and their nuclear luminosity evolution 

We define J\fL{ti, Mbh,o, L,t) {t > ti) so that J\fL{ti, Mbu,o, L,t)dtidMBH,odLdt is the comoving number 
density of local BHs with the following properties: the nuclear activity due to accretion onto their seed BHs^ 
was triggered during cosmic time ti ti+dti, the nuclei of their host galaxies were active and had luminosity 
in the range L ^ L + dL a,t cosmic time t, and the BHs are quiescent and have mass in the range Mbh,o ^ 
-^BH,o+dMBH,o at the present time to. Thus, the comoving number density of those BHs whose host galaxies 
have velocity dispersions in the range cr — > cr + do" is J\fL{ti: A'/bh.O; L, f)p((T|MBH.Oj to)diidMBH,odidtd(T [see 
p(c7|Mbh,0i ^o) in eq. 3]. We define the PDF of the nuclear luminosity of these BHs at cosmic time t whose 
host galaxies have velocity dispersion a as follows: 

^ Jodti J dMBH,o ^^L{t^,MB^ifi,L,t)p{a\MB^ifi,to) 

^ ~ JdL f*dti J dMBH,0 A'i(<„MBH,0,i,t)p((T|AfBH,0,io)' 

The function ^^{ti, Mbk,Oj L,t) is controlled by both the rate of triggering accretion onto seed BHs and 
the luminosity evolution of individual triggered nuclei. The rate of triggering nuclear activity is usually 
believed to be related to the formation and evolution of galaxies and is a function of cosmic time, which is 
very uncertain and not easy to predict. The luminosity evolution of individual triggered nuclei is believed 
to contain information on the accretion process in the vicinity of BHs and is a function of the physical time 
that the nuclei have spent since the triggering of the accretion onto seed BHs. Note that the luminosity 
evolution discussed here is different from the evolution of the characteristic luminosity of the QSO population 
as a function of redshift, which increases with increasing redshift at z < 2 — 3 and the variation tendency 
at 2; > 2 — 3 is not yet clear (e.g., see Fig. 6 in Boyle et al. 2000). The luminosity evolution is also not the 
evolution of the comoving number density of the QSO population as a function of redshift: the comoving 
number density of QSOs brighter than a certain luminosity has a peak at redshift 2; ~ 2 — 3 and decreases 
at both higher and lower redshift. As with AfL{ti, Mbb.ot L,t), the PDF p{L\a,t) is partly related to the 
triggering history and is not easy to predict. By integrating both the numerator and the denominator in 
equation (4) over the cosmic time t, we define a new nuclear luminosity PDF of the BHs with host galactic 



^Here the original mass of seed BHs does not come from gas accretion (which appears as QSO/AGN phenomena). Seed BHs 
could be remnants of Population HI stars, products of dynamical processes in dense star clusters, or primordial BHs formed 
in the early universe, etc. (e.g., van der Marcl 2004). The seed BH mass could also be due to non-luminous accretion. If the 
nuclear activity of a QSO/AGN is triggered by gas accretion recurrently, only the BH at the time of the first-time triggering 
(ti) is taken as the "seed BH" of the QSO/AGN. 
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velocity dispersion a as 

p.^l s ^ C /o dtj J dMsH.o ^Lft», Mbh.o, L, ^)p(a|MBH,o, ^o) 

~ / dL dt /o dii / dMBH.o ^^L{U, Mbh.o, t)p(a|MBH,o, to) ■ 

As is shown below, to obtain the PDF 'P{L\a), the necessary model inputs are only the luminosity evolution 
of individual triggered nuclei and the local Mbh.o-c relation [see also ^(MbhIct) similarly defined in eq. 23 
below] . 

As done in Yu & Lu (2004), for BHs with the same mass Mbh.o at present, we assume that their nuclear 
luminosity, >C(Mbh,o, ''") is a function only of the age of their nuclear activity r = t — ti. Then we use 
£(Mbh,0) t) to define two functions to describe the luminosity evolution of individual triggered nuclei. One 
is the lifetime of the nuclear activity for a BH with mass Mbh.o at present, defined by 



riife(MBH,o) = / dr (6) 

(7) 



. |d/:(MBH,0,T)/dT|, 

(see eqs. 13 and 14 in Yu & Lu 2004), where the integration in equation (6) is over the period that the nucleus 
was active and Tk(L, Mbh,o) {k = 1, 2, ...) in equation (7) are the roots of the equation £(Mbh,o, t) — L = 
(0 < T < to — U)- Given L and Mbh,Oi the number of the roots Tk can be more than 1, since £(Mbh,o,''') 
can be a non-monotonic function of r. The other function is a PDF of the nuclear luminosity, defined by 

" nue{Lu,o) ? |d£(MBH,o,r)/drk=.J 

so that Tiik{MBii,o)P{L\MBK,o)dL is the time that a BH (with mass Mbh,o at present) spent with nuclear 
luminosity in the range L ^ L + dL. According to equation (17) in Yu & Lu (2004), J\fL{ti, Mbka, L,t) is 
connected with the two functions defined above as follows: 

/ dt dtiJ\fL{U,MBH,0,L,t) = me{MBH,o)P{L\MBH,o)nMnH{MBH,0,to)- (9) 

Jo Jo 

Substituting equation (9) into equation (5), we have 

,0 ''■iifo(MBH,o)-P(i|A^BH,o)f^MBH(^BH,Oj to)p(o'|-MBH,Ojto) f^„-. 
J dAfBH.O Tiifc(AifBH,o)»^MBH(A'^BH,0,io)p(o-|MBH,0,io) 
/ dMBH.O Tiife(MBH,o)-P(^|MBH,o)p(MBH,ok,^o) 
J dMBH ,0 Tlife(-^BH,o)p(-^BH,o|o')io) 

where equation (3) is used. 



3.1.2. Observational counterparts 

The PDFs oip{L\(j. t) and 'P{L\a) can also be obtained directly from observations of QSOs/AGNs. Assuming 
that all the local massive BHs have experienced QSO/AGN phases, the observed QSOs/AGNs at rcdshift 
z may represent the progenitors of local galaxies and their central BHs at the cosmic time t{z), and the 
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QSO luminosity function (QSOLF) ^L{L,t) [defined so that ^LiL,t)dL is the comoving number density of 
QSOs/AGNs with luminosity in the range L ^ L + dL at cosmic time t] is given by (see also Yu & Lu 2004) 

^l{L, i) = ^ dt, J dMBH.O ^^Lit^, Mbh,0, L, t). (12) 

We can also define the QSO/AGN luminosity and velocity dispersion function (QSOLVF) $Ly(L, a, t) so that 
'^i,v{L, fT, t)dLda is the comoving number density of QSOs/AGNs with luminosity in the range L ^ L + dL 
at cosmic time t and galactic velocity dispersion in the range cr — > tr + dcr; and we have 

^Lv{L,a,t) = dti J dMBKfi ML{ti,MBH,o,L,t)p{a\MBH,o,to) (13) 

<ifL{L,t) = J <i>Lv{L,a,t)da. (14) 
Substituting equation (13) into equations (4) and (5), we have 

/-rl i\ ^Lv{L,(J,t) 



f^°dt<^>Lv{L,a,f) 



nm = /^l..? (16) 



3.1.3. V{a\L) 

Similar to the definition of 'P(Llcr), we may also define a velocity dispersion PDF in QSOs/AGNs at a given 
nuclear luminosity L as follows: 

1^. ^ Io° d* Jo dti / dMsH.o Ml{U, Mbh.o, L, t)p{a\MBH,o, h) 

JdaJ*''dtJ^dtiJdMBH,o^fLiti,MBn,o,L,t)pia\MBll,o,to) 

J dMBH,0 m^{MBHfi)PiL\MBH ,o)^Mbh (-^BH,Oj to)p(CT|MBH,0) to) ,^ „x 

/ dAfBH.O riife(MBH,o)i^(i'|MBH,o)nMBH(^BH,0,to) " 

SubstitTiting equations (13) and (14) into equation (18), we have the observational counterpart of V{a\L) 
given by 

/J" dt ^Lv{L,a, t) 
J^°dt^L{L,t) 

(see also eq. 64 in Yu & Lu 2004). 



As a summary of § 3.1, equation (16) or (19) gives the nuclear luminosity versus velocity dispersion rela- 
tion in QSOs/AGNs directly obtained from observations. Equation (11) or (18) shows the model prediction 
from the local Mbh.o-o" relation and the luminosity evolution of individual QSOs/AGNs. 



3.2. The BH mass versus velocity dispersion relation in QSOs/AGNs 



Given the luminosity evolution of an individual QSO/AGN £boi(-^BH,Ojt) (where the subscript "bol" repre- 
sents the bolometric luminosity) and the mass-to-energy conversion efiiciency e, the BH mass in QSOs/AGNs 
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follows the evolution below: 



MBH(MBH.o,r) = Mbh.(Mbh,o) + fdrii^lil^^MMMll), (20) 

Jo 



or 



A^BH(MBH,o,-r) = Mbh.o - / dr^ '- \ ■ (21) 

where e is assumed to be a constant and MBH.i is the seed BH mass. 

As done in § 3.1, we can also define the BH mass PDF by replacing L and £(Mbh,0)''") with Mbh and 

A^bh(Mbh,o, t). First we define AfMBuiky A'/bh,o, Mbh, (t > U) so that 

MMeniti, Mbh,0) MBH,i)dtidMBH,odMBHdf is the comoving number density of local BHs with the following 
properties: the nuclear activity due to accretion onto their seed BHs was triggered during cosmic time 
ti ^ ti + dti, the nuclei of their host galaxies were active with central BH mass in the range Mbh 
Mbh + dMBH at cosmic time t, and these BHs have mass in the range Mbh,o — > Mbh.o + dMBH,o at present 
time to. Similar to equations (4) and (5), we can then use A/mbh (^i> Mbh,o, Mbh, i) to define the BH mass 
PDF in the nuclear-active progenitors of local BHs, for example, 

I _ /odii/dMBH,0 AAMBH(^^,MBH,0,MBH,^)P(CT|MBH,0,^0) 
P[M-BK\0',t) = J (ZZj 

/ dMBH Jo dti J dMBH.O J^MuniU, MBKfi, MBli,t)pia\MBH,0,to) 



V{Mbm\(t) = /o° dt /o dti J dMBH.o TVmbh {ti, Mbh.o, Mbh, t)p(c7|MBH,o, ^o) ^^3) 

/ dMBH Jo° dt J* dti J dMBH.o J^M^uiU, Mbh.o, Mbh, i)p(CT|MBH,o, ^o) 

/ dMBH ,0 Tiifo(MBH,o)^'(MBH|MBH,o)p(MBH,o|c, to) 



J dMBH ,0 Tiife(MBH,o)p(MBH,o|f, to) 



(24) 



where 



dMRH 1 

P(MBH|MBH,o)dMBH = ,J rrr-, , (25) 

Tlife(MBH,o) |dA^BH(MBH, T)/dr|^^^/ 

and Ti is the only root of the equation A^bh(Mbh,o, t) — Mbh = (0 < r < to—ti) [note that A4bh(Mbh,o, t) 

in eq. 20 or 21 is a monotonically increasing function of r]. The BH growth history (or Tiifo) may vary for 
different Mbh,o- For simplicity, we assume that nife is independent of Mbh,o below, and thus equation (24) 
becomes 

V{MBYi\(y)= yP(MBH|MBH,o)p(MBH,ok,io)dMBH,o. (26) 

We can use moments to characterize the distribution of 7-*(logMBH|f ) = ln(10)MBH'P(MBH|cr), such as 
the mean of log Mbh at a given a, defined by 

(logMnnk) = j (logMBH)?'(logMBHk)dlogMBH, (27) 

and also the standard variance, the skewness, etc. The difference between the distributions or moments of 
■p(logMBH|f) and p(logMBH,o|o', io) contains information on BH accretion history and/or galaxy evolution. 
Obviously, by applying equation (26) to equation (27), we have the difference of the means of the distributions 
V{Mbb\o) and p(MBH,ok, to), 

^(logMBHk) = (log Mbh I (7) - (logMBH.ok) < (28) 
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[note that P(MBH|AfBH,o) = if Mbh > -^bh.o, and see (logMBH.ol'^) in eq. 1]. A positive difference of 
the means might suggest that galaxy velocity dispersions should increase during or after the nuclear activity 
(i.e., the assumption that the velocity dispersion does not significantly change during or after the nuclear 
activity should be revised). 

Similar to defining the QSOLVF as ^Lv{MBH,cr,t) in § 3.1, we can define the BH mass and velocity 
dispersion function in QSOs/AGNs (QSOMVF) as $mv(A'/bh, cr, t), so that $My (-^bh, cr, t)dMBH_<io- repre- 
sents the comoving number density of QSOs/AGNs with central BH mass in the range A^bh AfBH + dMsH 
and host galactic velocity dispersion in the range ct + da at cosmic time t. Thus we have 

"^""^"'"'^^^ /dMBH^MWMBH.a.t) ^''^ 
JdMBKJo°dt^Mv{MBK,C7,t) 

The distributions of p(Mbh|o', t) and V{Mbh\(^) give the BH mass versus velocity dispersion relations for 
QSOs/AGNs at a given cosmic time and for QSOs/AGNs at all redshifts, respectively. The V{Mbh\<7) is 
the Mbh-o" relation in QSOs/AGNs that is mainly discussed in this paper. 

Comparison of the expectation from equation (26) (or [24], [11], [18]) with observations of QSOs/AGNs 
(see eq. [30], [16], or [19]) can provide feedback to our understanding of the luminosity evolution of individual 
QSOs/AGNs, the mass growth of individual BHs, and the evolution of host galaxy velocity dispersions. 



4. The QSO/AGN luminosity evolution and BH growth model: £boi(AfBH,0) t) and 

A^bh(Mbh,o,t) 

In this section we present the detailed form of the luminosity evolution and BH growth [£boi(A^BH,0) ''') and 
A^BH(AfBH,0) t)] incorporated in the analysis in § 3. We assume that the luminosity evolution of individual 
QSOs/AGNs can be described by two phases (see also § 2.4 in Yu & Lu 2004): (1) after the accretion onto a 
seed BH is triggered, initially there is sufficient material to feed the BH, and the BH accretes with Eddington 
luminosity; (2) with the mass growth of the BHs and the consumption of the material to feed it, the material 
becomes insufficient to support Eddington accretion, and then the nuclear luminosity declines and is fainter 
than the Eddington luminosity (see also Yu & Lu 2004) . For simplicity, below we assume that the two phases 
appear only once for each BH (see also Yu & Lu 2004). 

We assume that the first phase lasts for a period of tj, and the BH mass increases to Mbh,i at time 
t = ti + Tl = tl. Thus, the nuclear luminosity in the first phase increases with time as 

£boi(r) = iEdd(MBH,i) exp ( ^^^^ ] < r < n, (31) 

where iEdd(-^BH,i) is the Eddington luminosity of a BH with mass Mbh.i, and 



Tsp = 4.5 X lO'^ 



0.1(1 -e) 



yr (32) 



is the Salpeter time (the time required for a BH radiating at Eddington luminosity to e-fold in mass). The 
BH mass increases as 

Mbu{t) = Mbh,i exp (^^] < r < n. (33) 



-lo- 



in the second phase we assume that the evolution of the nuclear luminosity declines as 
-CboKT) 



I LEdd(AfBH,i) exp (^-^^^ , for n < T < n + ^td, 

[O, forT>Tl+^TD, 



(34) 



where td is the characteristic declining timescale of the nuclear luminosity. We assume that QSOs/AGNs 
become quiescent when the nuclear luminosity declines by a factor of ?7 = exp(— ^) compared to the peak 
luminosity iEdd(-^BH.i), so there is a cutoff of the nuclear luminosity at t = tj + ^td in equation (34). The 
factor ^ is set to ln(10'^) = 6.9 (as in Yu & Lu 2004). According to equations (20) and (34), the BH mass 
increases as 



A^BH(r) = Mbh,i ^ 1 + — 



1 — exp 



Tl — T 



n<T<Tl+ ^Td. 



(35) 



With the assumption that the nuclear activity of all QSOs/AGNs is quenched at present (i.e., to — ti — ti ^ 
Td), the BH mass at present is given by 



Mbh.o = { 1 + ) Mbh.i. 



(36) 



Some constraints on £boi(-MBH,o, t) have been obtained by comparing the time integral of the QSOLF 
in observations with the prediction from the local BHMF in Yu & Lu (2004). For example, for the nuclear 
luminosity evolution models above, we should have n > rsp if e ^ 0.1 and td = and have rj > 0.2rsp if 
e ^ 0.31 and td = 0. The characteristic declining timescale of the second phase, td, should be significantly 
shorter than rsp, and BH growth should not be dominated by accretion in the second phase. In this 
paper, based on those constraints, we assume two family models for the parameter td: (1) td = and (2) 
Td = 0.3Tsp. In each of the family models, we consider four cases for the total mass increase of a BH during 
its nuclear activity period: that is, Mbh.o / Mbka = exp(0.3), exp(l), exp(2), and exp(4), which are denoted 
as cases A-D, respectively. Thus, we have td = and ri/rsp = 0.3, 1, 2, and 4 for models lA-lD; and we 
have td/tsp = 0.3 and n/rsp = 0.04, 0.7, 1.7, and 3.7 for models 2A-2D. 



5. Simulation results of the Mbh-c relation in QSOs/AGNs 

In this section, by applying the BH growth models in § 4 to the analysis in § 3, we obtain the expected 
observational appearance of the Mbh-c relation in QSOs/AGNs. We quantitatively illustrate how this 
relation in QSOs/AGNs depends on the BH growth history and how the flux limit of telescopes affects 
the observational appearance of the relation. Similarly, the expected observational appearance of the L-a 
relation in QSOs/AGNs can also be obtained, but for simplicity we do not show and discuss it in this paper. 



5.1. n and td 

With the local Mbh.o-ct relation [or p(Mbh,o|o', ^o); see eqs. 1 and 2] and the QSO/AGN luminosity evo- 
lution and BH growth models described in § 4, we use equation (26) to get the distribution V{Mbh\<^) in 
QSOs/AGNs. The results of P(logMBH|o' = 200 kms^-'^) arc shown in Figure la for model 1 and in Figure lb 
for model 2. In each panel we show the distribution in nearby normal galaxies, p(log Mbh,o|C) io)) by the 
solid line and show the distribution in QSOs/AGNs, 'P(logMBH|o'), by the dotted, short-dashed, dot-dashed, 
and long-dashed lines for cases A-D, respectively. Note that the shape of the distribution p(logMBH,o|o', fo) 
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has been assumed to be independent of a, and the characteristic increasing and decreasing timescales (ti, 
td) have been assumed to be independent of Mbh,o; therefore, the shapes of the curves of P(logMBH|o') are 
independent of cr, and the distributions for other velocity dispersions can be obtained simply by shifting the 
curves to higher or lower BH masses. The difference of the distributions ■p(logMBHlCT) and p(logMBH,o|o', to) 
can be characterized by the difference of the mean (5(log Mbh|c) in eq. 28), the standard variance, the skew- 
ness, etc. For simplicity, wc mainly discuss the difference of the mean below. 

As seen from Figure 1, with increasing timescale tj (see eq. 31), the curve peak of 7'(logMBH|'7) shifts to 
low BH mass and the scatter of ^(log MbhI^) increases, since the probability of observing QSOs/AGNs with 
low BH masses increases. Our calculations show that we have the difference 5 (log Mbh) = log(MBH,i/AfBH,o)/2 — 
-0.07,-0.2,-0.4 and ,-0.9 dcxfor models lA-lD in panel (a) and (5(logMBH) ^ -0.02,-0.08,-0.2, and 
—0.6 dexfor models 2A-2D in Figure lb. The difference (5(logMBH) is difficult to detect for cases (A) and 
(B), considering the current large measurement errors of BH masses in QSOs/AGNs (see the discussion in 
§ 6), and the difference is relatively significant for case D, which should be detectable if there existed a 
sufficiently large and complete sample of QSOs/AGNs. As seen from Figure 1, for the same case of the 
two family models, compared to model 1 (td = 0), 7^(logMBH|cr) m model 2 (td = O.Srsp) is lower at the 
low-mass end and higher at the high-mass end, and the peak of the curve is closer to the peak of the local 
distribution p(log Mbh.oIc, ^o), although the BH mass is increased by the same factor during the nuclear 
activity. The reason is that in model 2, QSOs/AGNs have a substantially large probability of being observed 
in the second phase (see eq. 34), in which the QSOs/AGNs are shining at sub-Eddington luminosities and 
BH masses grow relatively slowly and have acquired nearly all of their mass. In addition, Figure 1(b) shows 
that the existence of sub-Eddington accretion at the late evolutionary stage may also introduce skewness to 
the distribution of ^(logMBHlf ). 



5.2. Super-Eddington accretion 

In Figure 1 the nuclear luminosity has been assumed to be limited by the Eddington luminosity. In the 

literature, however, it is argued that QSOs/AGNs might be able to accrete material at a rate higher than 
the Eddington accretion rate (which is conventionally defined by the Eddington luminosity for e = 0.1 
here) at the early stage of their evolution (since initially the material supply for the BH growth may be 
sufficiently large; e.g., Blandford 2004). Below we use a simple example to illustrate how the Mbh-c relation 
in QSOs/AGNs is affected by super-Eddington accretion. 

We modify model ID (ri/rsp = 4) as follows: BHs first accrete material at a rate higher than the 
Eddington rate by a factor of I in the initial 3Tgp period of the nuclear activity, where I is a constant and 
Tgp = l~^TSp (eq. 32), and then BHs accrete material at the Eddington rate for a period of rsp. We denote 
the modified model as model ID'. The BH mass is increased by the same factor in the initial Srgp period 
of the nuclear activity in model ID' as that in the initial 3Tsp period of the nuclear activity in model ID. 
The total BH mass increase over the whole nuclear activity period is also the same in models ID' and ID. 
The ■p(logMBH|o') of model ID' are shown as the dotted {I = 2) and dot-dashed {I = 10) lines in Figure 2. 
For comparison, p(logMBH,o|f, to) and the results of models IB (ri/rsp = 1) and ID shown in Figure la are 
also shown in Figure 2 as the solid, short-dashed, and long-dashed lines, respectively. As seen from Figure 2, 
the dot-dashed line (model ID' with I = 10) is lower than both the dotted (model ID' with / = 2) and 
long-dashed lines (model ID) at the low-mass end (< lO'^M©), and the distribution of the dot-dashed line is 
closer to that of the short-dashed line (model IB; tj/ rsp = 1) than that of the dotted line. The reason is that 
the larger the accretion rate, the more rapidly the BH mass increases, and thus the smaller the probability of 
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the low-mass BHs within a unit mass range being observed, since the high-rate or super-Eddington accretion 
is at the early BH growth stage here. For example, in model ID', the BH mass increases by a factor of ~ 20 
at the early super-Eddington accretion stage and by a factor of --^ 2.7 at the Eddington accretion stage; for 
the case of Z = 2, the time fractions for a QSO to be at those stages are 60% and 40%, respectively, but 
for the case of / = 10, the corresponding time fractions become only 23% and up to 77%. The effect of the 
super-Eddington accretion can also be quantitatively characterized by high-order moments of P(logMBH|o'), 
such as skewness. 



5.3. Effects of the flux limit of telescopes 

In the results above, we have assumed that QSOs/AGNs with any luminosities can be observed, without 
being limited by the flux limit of telescopes. In this subsection we show how the observational appearance 
of the Mbh-o" relation in QSOs/AGNs is affected if only those QSOs/AGNs brighter than a certain absolute 
magnitude are detected. We show the results of ■p'(logMBH|o'), the portion of 'P(logMBH|o') contributed 
to QSOs/AGNs with such absolute magnitude truncations, in Figure 3. Unlike those of V {\og MbhI^') in 
Figures 1 and 2, the shapes of 7-"(logMBH|cr) vary for different galactic velocity dispersions because of the 
magnitude truncation. Thus, we show the results at two different velocity dispersions: a = 200kms~^ in 
Figure 3a and 3b, and a — 320kms~^ in Figure 3c and 3d. For a = 200kms~^ we set the magnitude trun- 
cation of the QSOs/AGNs to be Mb < —20 and —23 (the bolomctric correction for the B-band luminosity, 
Cb, is defined through Lboi = CbLub^ where is the energy radiated at the central frequency of the B 
band per unit time and logarithmic interval of frequency, and we set Cb = 11.8; see Elvis et al. 1994). For 
the high-velocity dispersion a = 320kms~^, we set a brighter magnitude truncation {Mb < —23 and —25) 
because of the following considerations: at z < 2-3 the characteristic luminosity of QSOs/AGNs increases 
with the increase of redshift (see the QSOLF given in Boyle et al. 2000), or high- mass BHs are more likely 
to grow up at high redshifts, and for high-redshift objects, only the bright ones are observed. Other model 
parameters in Figure 3 (tj, td) arc the same as those in models 2A-2D (sec Fig. lb). As seen from Figure 3, 
QSOs/AGNs housing small BHs cannot be detected because of the faintness cutoff. The distribution of 
■p'(logMBH|c) is also lower than that of V{\ogMBii\a) at the high- mass end because some QSOs/AGNs 
that house big BHs but are at the late stage of the second luminosity evolutionary phase (eq. 34) may 
also be missing from observations. Our calculations show that if the sample of QSOs with Mb < —20 is 
complete, |(5(logMBH)| is about 0.6 dexat a = 200kms~^ for model 2D (see the long-dashed line in Fig. 3a), 
which is detectable if the accuracy of the BH mass measurement is within a factor of ~ 2-3. However, if the 
magnitude truncation is Mb < —23, we have the difference |(5(logMBH)| ^ 0.1 dex at a — 200 km s"""^ (see 
Fig. 3b), which is too small to detect currently. In short. Figure 3 suggests that the flux limit of telescopes 
may weaken the difference of the observed Mbh-c relation in QSOs/AGNs from the local Mbh,o-o' relation, 
since small BHs, the main contributors to the difference of these relations, may be missing from observations. 

To further illustrate the appearance of the MBH-f relation in QSOs/AGNs, we simulate a sample of local 
galaxies and QSOs/AGNs and show the result in an Mbh versus a plot (see the black dots for local galaxies 
and red dots for QSOs/AGNs in Fig. 4). To do this, we flrst use the Monte-Carlo method to draw a sample of 
10^ galaxies from the velocity dispersion distribution of local galaxies with a > 80 kms^^ (see the black dots 
in Fig. 4; for the velocity dispersion distribution of local galaxies, see Yu & Lu 2004). Then, for each (local) 
galaxy, we use the Monte-Carlo method to generate one progenitor (i.e., QSOs/AGNs) from the distribution 
function P{MBii\a). Note that the number of progenitors generated for each local galaxy, which should be 
proportional to the lifetime of the nuclear activity of this galaxy rufe, is the same, since the lifetime has 
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been assumed to be the same for all the galaxies here. In addition, although only one progenitor is produced 
for each galaxy, the total number of progenitors generated for all the galaxies is large enough in this paper 
to illustrate the appearance of the Mqu-'ict relation in QSOs/AGNs. Finally, we cut off faint progenitors, 
and the rest are just the sample of QSOs/AGNs with magnitude truncation (red dots in Fig. 4). We show 
the results of models 2C and 2D (td/tsp = 0.3, Ti/rgp = 2,4) with an absolute magnitude truncation of 
Mb < —20 in Figures 4a and 4b, and show the results of models 2C and 2D with Mb < —23 in Figures 4c 
and 4d, respectively. The simulation results of models 2A and 2B are not shown in this figure because they 
are visually difficult to distinguish from the local Mbh.o-ct relation. In models 2C and 2D, the increase of 
the BH mass via accretion is large enough (e.g., by a factor of > 10) that as seen from Figures 4a and 4b, if 
all or at least most of the QSOs/AGNs with absolute magnitude Mb < —20 can be detected, many red dots 
with small BH masses (simulated QSOs/AGNs) will be below the black dots (simulated local galaxies), and 
the M^ii-a relation in QSOs/AGNs will be easily distinguished from the local M^y^ Q-a relation. However, 
the Mbh ^ relations with truncation Mb < — 23 in Figures 4c and 4d are not as easy to distinguish from 
the local Mbh.o-c relation as those in Figures 4a and 4b, since QSOs/AGNs housing small BHs are missing 
from observations. 

In addition, the observational appearance of the Mbh-o" relation in QSOs/AGNs may also be affected 
by obscuration. If QSOs/AGNs at the early evolutionary stage are strongly obscured, as suggested by 
Fabian (1999) and King (2003), the Mbh-c relation in the observed QSOs/AGNs should be less deviated 
from the MBH,o-f relation in nearby galaxies, since the early stage of BH growth is missing from obser- 
vations. However, if the obscuration is a purely geometric effect (as suggested by the unification model of 
QSOs/AGNs) and the obscured fraction does not depend on the QSO/AGN luminosity or central BH mass, 
the observational appearance of the Mbh-o" relation will not be affected by obscuration. 

6. Discussion on current observational results of QSOs/AGNs 

Some observations of the BH mass/nuclear luminosity versus velocity dispersion relation in QSOs/AGNs 
have been made in the past several years (e.g., Laor 1998; Wandel 1999; Gebhardt et al. 2000b; Ferrarese ct 
al. 2001; Shields ct al. 2003). For example, Gebhardt et al. (2000b) and Ferrarese ct al. (2001) argue that 
QSOs/AGNs follow the same BH mass versus velocity dispersion relation as local galaxies based on some 
galaxy samples including a small number (seven or six) of QSOs/AGNs that have measurements of both 
BH masses and velocity dispersions. The masses of BHs in QSOs/AGNs in these studies arc estimated by 
the reverberation mapping technique (e.g., Wandel et al. 1999; Kaspi et al. 2000). However, Krolik (2001) 
argues that the BH mass measured by this method may be either underestimated or overestimated by a 
systematic error of up to a factor of 3, partly because of the deviation of the simple assumptions on 
the dynamics/kinematics and the geometry of broad- line regions from the reality. For another example, 
Shields et al. (2003) use a larger sample of QSOs/AGNs (~ 100, including ~ 20 high-redshift [z > 1] QSOs) 
with estimated BH mass and velocity dispersion and also argue that the BH mass versus velocity dispersion 
relation in QSOs/AGNs is consistent with the relation in nearby galaxies. In Shields et al. (2003), the BH 
mass is estimated through the empirical law between BH mass and broad-line region size (e.g., Kaspi et 
al. 2000) and the [OIII] line width is used as a surrogate for the galactic velocity dispersion (e.g.. Nelson 
2000). These methods might introduce further uncertainties in BH mass and velocity dispersion. Despite 
the measurement error (or possible systematic error) in BH mass and velocity dispersions, we note that the 
flux limit of telescopes can still cause consistency between the relations in Shields et al. (2003), especially 
for the subsample of high-redshift QSOs (which include only BHs with masses higher than 10^ M©). Just 
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log(MBH/Me) log(MBH/M^) 

Fig. 1. — Estimated BH mass distribution in QSOs/AGNs [i.e., ■p(logMBH|o'); see eq. 26] at a given host 

galaxy velocity dispersion a ~ 200 km s""'^. The luminosity evolution and BH growth models used are 
described in § 4. In (a) are the results for model 1 (td = 0) and (b) is for model 2 (ro/rsp = 0.3). In 
both panels, the solid line represents p(logMBH,o|o', to) (or the local MBH,o-f relation, for which a Gaussian 
distribution of logMeH is assumed with an intrinsic scatter of 0.27 dex; see eqs. 1 and 2). The dotted, 
short-dashed, dot-dashed, and long-dashed lines represent the results obtained by setting ln(MBH.o/AfBH.i) = 
0.3, 1, 2, and 4 (which are denoted by cases A-D, respectively). That is, we have ri/rsp = 0.3, 1, 2, and 4 for 
models lA-lD in (a), and have ri/rsp = 0.04,0.7, 1.7, and 3.7 for models 2A-2D in (b). As seen from this 
figure, with increasing t\ the probability of ^(logMBHlCT) at the low-mass end increases, and the deviation of 
■p(logMBH|o') from p(log Mbh,o|c, to) increases. Compared to the results of the purely Eddington accretion 
model in (a), the existence of sub-Eddington accretion at the late evolution stage (td 7^ in [b]) decreases 
the difference |(5(logMBH)| (eq. 28) and also may introduce skewness to the distribution of ■p(logMBH|o'). 
See details in § 5.1. 
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log(MB„/Mj 

Fig. 2. — Estimated BH mass distribution V {log MbhIct) in QSOs/AGNs at a given host galaxy velocity 
dispersion, for which a super-Eddington accretion at the early evolutionary stage of the nuclear activity is 
assumed. The solid, short-dashed, and long-dashed linos arc the same as those in Figure 1(a). The dotted 
and dot-dashed lines represent the results obtained by model ID', assuming that BHs accrete material at 
the Eddington rate only at the final rsp period of the nuclear activity. Before that period, BHs accrete at a 
rate higher than the Eddington rate, for example, by a factor of 2 {dotted line) or 10 {dot-dashed line). As 
seen from this figure, supcr-Eddington accretion at the early stage decreases V {A'Ibh\o') at the low-mass end 
and may also introduce skewness to the distribution V (log Mbh\<7) ■ See details in § 5.2. 
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Fig. 3. — Estimated (log Mbh I cr), i.e., the portion of P (log MbhIo") contributed to by QSOs/AGNs brighter 
than a certain absolute magnitude. The line types have the same meanings as in Figure 1. The i?-band 
absolute magnitude truncations are Mb < —20,-23, and —25 in (a), (b-c), and (d), respectively. Other 
model parameters (ti, td) are the same as in model 2 (see Fig. lb). Panels a and b are for the distributions 
at fj = 200kms~^, and panels c and d are for the distributions at tr = 320kms~^. This figure (and also 
Fig. 4 below) suggests that the flux limit of telescopes may decrease the difference |(5(logMBH)| between the 
observed Mbh-c relation in QSOs/AGNs and the Msufi-o' relation in local galaxies, since small BHs, the 
main contributors to the difference, may be missing from observations. See details in § 5.3. 
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Fig. 4. — Simulation results of the Mbh-c relation in QSOs/AGNs. The black dots represent local galaxies 
with an assumed intrinsic scatter in logMeH.o of 0.27 dex, and the red dots represent QSOs/AGNs. Panels 
a and c are for model 2C, and panels b and d are for model 2D. Panels a and b only include QSOs/AGNs 
with AIb < —20, and panels c and d only include QSOs/AGNs with Mb < —23. Each panel has 10^ black 
dots; and the numbers of red dots are 4276, 3562, 658, and 462 in panels (a-d), respectively. See details in 
§ 5.3. 
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as illustrated in Figures 3 and 4, since low-mass progenitors of big BHs are not detected, the relation in 
QSOs/AGNs is difficult to distinguish from that in local galaxies. 

Ignoring possible (systematic) measurement errors in BH masses and velocity dispersions or selection 
effects for the samples of QSOs/AGNs in Gebhardt et al. (2000b) and Ferrarese et al. (2001), we give a 
tentative discussion below on the possible implications of their results. The BH masses in QSOs/AGNs have 
a negative offset of —0.21 dex from local BHs for the sample in Gebhardt et al. (2000b) (including seven 
QSOs/AGNs) but have a positive offset for the sample in Ferrarese et al. (2001) (including six QSOs/AGNs). 
According to equation (28), a positive offset might mean that the velocity dispersion increases during or after 
the nuclear activity. However, here we think that the difference of the offsets for these two different samples 
is partly because (1) the slope of the local Mbh.o-c relation (see 7 in eq. 1) in Ferrarese et al. (2001) is 
steeper than that in Gebhardt et al. (2000b), and (2) the velocity dispersions of two AGNs (NGC 4051 
and NGC 4151, which are included in both samples) measured in Ferrarese et al. (2001) are smaller than 
those quoted in Gebhardt et al. (2000b) by ~ 10%. According to Gebhardt et al. (2000b), the offset of the 
-^BH — o' relation in QSOs/AGNs from the Mbh.o — o' relation in nearby galaxies can be down to — 0.21dex. 
The amount of the offset, based on the calculation results illustrated in Figure 1(a) and (b), suggests that 
the period in which QSOs/AGNs accrete material with Eddington luminosity can be as long as ~ 1-2tsp and 
the BH mass can be increased by a factor of up to ~3-8 [i.e., ~ exp(l)-exp(2)] during the nuclear activity 
(this suggestion does not contradict the result obtained in Yu & Tremaine 2002 and Aller & Richstone 2002 
that local BH mass comes mainly from accretion during QSO phases). 

Since current samples of QSOs/AGNs with BH mass and velocity dispersion measurements suffer from 
the uncertainty of small number statistics, large measurement errors, and incompleteness due to the flux 
limit of telescopes, a large complete sample of QSOs/AGNs with accurate BH mass and velocity dispersion 
measurements is needed to obtain quantitative constraints on the BH accretion history from the comparison 
of the BH mass versus velocity dispersion relation in QSOs/AGNs with that in local galaxies. 

7. Conclusions 

In this paper we have investigated how the BH mass/nuclear luminosity versus velocity dispersion relation 

in QSOs/AGNs is connected with the local BH mass versus velocity dispersion relation through the mass 
growth/nuclear luminosity evolution of individual BHs (see eqs. 24 and 11), by ignoring BH mergers and 
assuming that the velocity dispersion does not significantly change during and after the nuclear activity phase. 
Circumventing the uncertain history of triggering the accretion onto seed BHs. the rcilation in QSOs/AGNs 
considered in this paper is for a sample of QSOs/AGNs at all redshifts, rather than for QSOs/AGNs at a 
given cosmic time. 

Using the observed local BH mass versus velocity dispersion relation and the assumed form of the 
QSO/AGN luminosity evolution and BH growth, we simulate the observational appearance of the BH mass 
versus velocity dispersion relation in QSOs/AGNs. The simulation results quantitatively show how the 
BH accretion history (e.g., characterized by the timescales tj and td in this paper; see § 4) affects the 
difference between the relation in QSOs/AGNs and that in local galaxies. A simple example to illustrate 
this is that: if the BH mass increases by a factor of > 10, mainly via Eddington accretion, the relation 
in QSOs/AGNs will significantly deviate from the relation in nearby galaxies, with the mean logarithm 
of BH mass at a given velocity dispersion being smaller than that of local BHs by an absolute difference 
of > 0.3 dex(see Figs, la and lb). We quantitatively show that if QSOs/AGNs are not always accreting 
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material at the Eddington rate (e.g., accreting at super-Eddington rates at their early evolutionary stages 
or at sub-Eddington rates at their late evolutionary stages), the distribution function of BH mass at a given 
velocity dispersion in QSOs/AGNs can be skewed (i.e., the probability of observing QSOs/AGNs with low 
BH mass will be relatively decreased compared to that for QSOs/AGNs with high BH mass). We also show 
that the observational difference between the relation in QSOs/AGNs and that in local galaxies may be 
weakened by the selection effect due to the flux limit of telescopes, since QSOs/AGNs housing small BHs 
may be missing from observations. To further constrain the BH accretion history from the BH mass versus 
velocity dispersion relation in QSOs/AGNs, a large, complete QSO/AGN sample with accurate BH mass 
and velocity dispersion measurements from observations is needed. 

In addition, the observational appearance of the nuclear luminosity versus velocity dispersion relation 
in QSOs/AGNs can also be predicted by the method described in this paper, and it would also be useful to 
compare the expectation with future observations. 
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